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The impact of superthermal electrons on dispersion properties of isotropic plasmas and on the 
modulational instability of a monochromatic Langmuir wave is studied for the case when the power- 
law tail of the electron distribution function extends to relativistic velocities and contains most 
of the plasma kinetic energy. Such an energetic tail of electrons is shown to increase the thermal 
correction to the Langmuir wave frequency, which is equivalent to the increase of the effective electron 
temperature in the fluid approach, and has almost no impact on the dispersion of ion-acoustic waves, 
in which the role of temperature is played by the thermal spread of low-energy core electrons. It 
is also found that the spectrum of modulational instability in the non-maxwellian plasma narrows 
significantly, as compared to the equilibrium case, without change of the maximum growth rate and 
the corresponding wave number. 
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I. INTRODUCTION 

The typical feature of turbulence evolution in a beam- 
plasma system is the formation of superthermal tails on 
the distribution function of plasma electrons. These tails 
are often observed in space and laboratory plasmas. In 
particular, the electron distribution in the solar wind (see 
the review paper pQ and references therein) was found to 
be well approximated by the family of kappa-functions 
[2HB], decaying in velocity as a power law. Moreover, 
it has been shown recently that these functions describe 
asymptotically stationary solutions of weak turbulence 
equations, in which self-consistent evolution of both wave 
and particles spectra is taken into account [7HS] • In labo- 
ratory experiments, beam-plasma interaction appears to 
be much more intense than in the solar wind, that is why 
it results in the formation of strongly non-maxwellian 
electron distribution, in which most of the plasma ki- 
netic energy is concentrated in a rather small population 
of fast electrons. Such a slowly decaying momentum dis- 
tribution of plasma electrons (/ cx p~ 5 ) was really ob- 
served in the experiments on turbulent plasma heating 
by the powerful electron beam in the multi-mirror trap 
GOL-3 [10]. 

It is obvious that such an energetic tail of superther- 
mal electrons should modify not only the linear disper- 
sion of plasma modes, but also probabilities of various 
nonlinear processes responsible for formation of the tur- 
bulent spectrum. In the experiments of interest [llj , the 
high-current relativistic electron beam interacts with the 
plasma in the regime of strong plasma turbulence. In this 
case, the main nonlinear process responsible for trans- 
ferring the wave energy of beam-driven plasma oscilla- 
tions to the nonrcsonant part of turbulent spectrum is the 
modulational instability. The goal of this paper is to in- 
vestigate how the dispersion of linear plasma modes and 
the growth rate of the modulational instability pumped 
by the monochromatic Langmuir wave are modified in 



a strongly non-maxwellian plasma. Our interest to this 
problem is motivated by the need to interpret the experi- 
mental information about the intensity and the frequency 
spectrum of electromagnetic radiation generated in the 
turbulent plasma near the second harmonic of plasma 
frequency. Recent calculations of the radiation power 
[L2l [13] for these experiments were based on the well 
known analytical model of strong plasma turbulence [14] . 
According to this model, the typical wavenumber of mod- 
ulational instability determines the size of the energy- 
containing region of turbulent spectrum, and the typical 
growth rate allows to estimate the saturation level of tur- 
bulence energy. Thus, the study of the modulational in- 
stability in a non-maxwellian plasma will clarify whether 
superthermal electrons are able to modify the main pa- 
rameters of emitting part of the turbulent spectrum. 

The linear analysis of Langmuir waves in an isotropic 
plasma with a power-law momentum distribution in the 
framework of relativistic kinetic theory has been applied 
recently to the solar wind plasma near the orbit of the 
Earth [15] . In our paper, we study the dispersion prop- 
erties of much hotter plasma that is typical to the lab- 
oratory beam-plasma experiments, and we focus on the 
validity of various widely used approximations. Inves- 
tigations of modulational instability in relativistic plas- 
mas with the Maxwell- Jiittner distribution jTHl [17] or 
in the kappa-distributed plasmas with nonrelativistic su- 
perthermal tails [THIHH] are a l so °f interest in last years. 
In contrast to these papers, we consider the modula- 
tional instability of a Langmuir wave in an unmagnetized 
plasma with the relativistic power-law momentum distri- 
bution of plasma electrons. 

In Section II we obtain numerical solutions of the lin- 
ear dispersion equation for the potential plasma waves 
of isotropic non-maxwellian plasma and compare them 
to the results predicted by fluid and kinetic approxima- 
tions. In Section III, for the modulational instability, 
we deduce the dispersion equation that takes into ac- 
count relativistic and kinetic effects for high-frequency 
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oscillations, and compare numerical solutions for unsta- 
ble spectra with the case of maxwellian plasma. In the 
concluding Section IV, we formulate our main results. 



II. LINEAR DISPERSION OF POTENTIAL 
PLASMA OSCILLATIONS 

Let us analyze the proper waves of non-maxwellian 
plasma assuming cold ions and using the following dis- 
tribution function for electrons: 



the whole distribution is 



/(P) 



Co H( Ph -p) 
47T (p 2 + Ap 2 ) 5 / 2 



(1) 



where H(p/ l — p) is the Heaviside step function, p^ is 
the threshold momentum, above which there are no elec- 
trons in the plasma, Ap is the typical momentum spread 
of plasma electrons, and Cq = 3Ap 2 / sm 3 (ph/Ap) is the 
coefficient corresponding to the normalization J fdp = 1. 
For the parameters Ap = 0.066m e c and ph — 5m e c (m e 
is the rest mass of electons, c is the speed of light), this 
function is best suited for the description of real electron 
distribution observed in beam-plasma experiments at the 
GOL-3 facility. The feature of this distribution function 
is that, for small momenta p < p = 2Ap, it does not differ 
greatly from the maxwellian function with the tempera- 
ture T = 1 keV and, for large momenta p > p, it has 
a power-law tail decreasing so slowly that most of the 
plasma kinetic energy is concentrated in superthermal 
electrons (fig. [l}. The tail component of this particular 
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Figure 1: Momentum distribution of plasma electrons /(p) 
for parameters typical to beam-plasma experiments at the 
GOL-3 facility. 

distribution (p > p) contains 28% of the particles and 
89% of the kinetic energy; the effective temperature of 



T e ff = ~f m e c 2 ( 7 - l)/dp = 6.5 keV, (2) 

whereas the temperature of core electrons with the rela- 
tive density n c = 72% turns out to be much smaller: 

J m e c 2 (7-l)/dp = 0.97keV (3) 



2 

3n c 



(7 is the relativistic factor). 

It is well known that dispersion laws of Langmuir and 
ion-acoustic waves in the maxwellian plasma have the 
following dimcnsionless forms: 



uj £ ~ 1 + -k 2 T, 



kT 1 ' 2 




(4) 



(5) 



Here frequencies are expressed in units of the plasma 
frequency uj p — (47re 2 n/m e ) 1/ ' 2 , where e is the electron 
charge and n is the plasma density, wavenumbers are 
measured in units of u> p /c, and the electron temperature 
in units of m e c 2 . These formulas have a simple hydro- 
dynamic interpretation, according to which the electron 
pressure gradient V^p = 7 e TV<5n drives low-frequency 
acoustic oscillations of ions ( 7e = 1) and creates an ad- 
ditional force for high-frequency Langmuir oscillations of 
electrons ( 7e = 3). To adapt these insights to the case 
of non-maxwellian plasma, we should clarify what tem- 
perature {T e ff or T c ) is more appropriate for the correct 
description of wave dispersion. 

As we shall show later, the same expressions for real 
parts of wave frequencies can be derived from the exact 
kinetic theory, if we use the hydrodynamic approxima- 
tion for Langmuir waves and kinetic approximation for 
ion-acoustic waves. From the nonrelativistic theory one 
can conclude that, for Langmuir oscillations, the role of 
temperature is played by 



T f = 



v 2 fdp = 



(6) 



whereas, for ion-acoustic oscillations, the temperature is 
defined in another way 



T s = 



1 



< 1/v 2 > ' 



(7) 



In the maxwellian plasma with the temperature T, these 
values are equal, Tg — T s = T, but even from these 
simple formulas one can see that relative contributions 
of superthermal electrons to Tg and T s differ substan- 
tially. In a non-maxwellian plasma, the energetic tail 
component should lead to the significant increase of Tg, 
since it is proportional to the averaged kinetic energy of 



3 



particles, and should not contribute to the value T s , in 
which the inverse energy should be averaged. Thus, in 
the nonrelativistic case, superthermal electrons increase 
the thermal correction to the Langmuir wave frequency 
that is determined by the effective temperature of the 
whole distribution T e ff and have almost no impact on 
ion-acoustic waves, for which the thermal spread of core 
electrons T c is more suitable to the role of temperature. 

Let us now find out how these results change, if the tail 
of superthermal electrons extends to relativistic energies. 
First we will derive dispersion laws w r (fc) and damping 
rates T(k) of potential plasma waves in the hydrodynamic 
and kinetic limits, and then investigate applicability of 
these approximations by comparing approximate results 
with numerical solutions of the exact dispersion relation. 

For slowly damping waves of isotropic plasma, the dis- 
persion equation 



1 



1 

fc 2 



k • df/dp 



kv 



dp 



m e /m t 







(8) 



can be reduced to the form 



Re£||(w r — ir, k) + ilm£ii(cj r , k) = 0, 

where real and imaginary parts of dielectric permittivity 
are defined by the integrals 

Rce \\ = 1 + T2 VP r - d P 2~ ' ( 9 ) 

k z I v u — kv lo z 



Imeii = ~ ( —^fs(cj - kv)rfp. (10) 

11 k z J v op 

The dispersion law in this case is found from the equation 
Re£||(cj r ,fc) =0, (11) 
and the damping rate is calculated as 
( I m£ || 



r 



\9Ree||/0a 



(12) 



We can calculate the electron contribution to Reey using 
the hydrodynamic approximation (u> 3> kv) for Lang- 
muir waves and the kinetic approximation (cj -C kv) for 
ion-acoustic waves. It means that, in the former case, we 
can use the expansion 



1 = 1 A , kv (kv) 2 (kvf 
lu — kv 



U) \ u> 
and in the latter case 

1 



to — kv 



1 

k^ 



(13) 



(14) 



Im Eii is determined for both types of waves by the unified 
expression: 



Imeii 



2ir 2 



k 3 



lofiPo) 



pfdp 



IHl 



(15) 



where po is the momentum of electrons getting into 
Cherenkov resonance with the wave, 



oj r /k 



and 70 = \f \ + p\ is the corresponding relativistic factor. 
Thus, the real frequency of Langmuir wave is governed 
by the equation 



k*T, 



where 



= 0, 



l-v 2 /3 



T.-t± l-U 



37 V 5 



(16) 



(17) 



(18) 



and the real frequency of ion-acoustic wave is determined 
by 



where 



1 m„ Irrii 
1 H — — - = 



rp — 1 / 



pv 



(19) 



(20) 



In the nonrelativistic limit, ujq — > 1 and calculation of 
electron temperature really reduces to averaging either 
the kinetic energy Tg — >< w 2 /3 > or the reciprocal value 
T^ 1 —>< 1/v 2 >. From these formulas, it is seen that 
relativistic effects manifest themselves in decreasing the 
plasma frequency (lJq < 1) due to weighting of tail elec- 
trons. Moreover, the temperature Tg in the relativis- 
tic case ceases to be proportional to the total kinetic 
energy and can no longer be identified with the value 
T e ff —< 2(7 — l)/3 >. Indeed, the expression averaged 
in Tg ( 18 ) differs from the value 2(7 — l)/3 by the factor 



F 



(7 + 1) 
57 3 



3 



(21) 



which substantially reduces the relative contribution of 
fast electrons to the thermal correction to the Langmuir 
wave frequency. As to the temperature T s , its calculation 
for both low 7 <C 1 and high 7> 1 energies reduces to av- 
eraging the value that is inversely proportional to the en- 
ergy. It means that relativistic effects do not change the 
earlier conclusion that superthermal electrons do not con- 
tribute to the electron pressure perturbation driving low- 
frequency acoustic oscillations. For the distribution func- 
tion ( T]) , the temperature Tg accounting for the relativis- 
tic effects appears to be twice smaller than the effective 
temperature of the whole distribution (Tg — 0.56 T e //) 
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and the temperature T s remains close to the tempera- 
ture of low-energy core electrons (T s — 1.15 T c ) as in the 
nonrelativistic case. 

The exact dispersion relation for the damping plasma 
modes can be written in the form 



e\\(w,k) = 1 



m e /m l 4tt f / 2 2 

-(7 +P )dp 



k 2 J 7 
o 



df ( u) — kv 



o 

where the imaginary part of dielectric permittivity is still 
defined by (15) and the function ln(...) means the prin- 



cipal branch of complex logarithm. For Langmuir waves, 
the comparison of numerical solution of this equation 
ui = ui r — iT with the results of the hydrodynamic ap- 
proximation 



/ koj t Cq 



1/2 



(23) 



6fc 3 (wg + 6k 2 T e /uj 2 ) (p§ + Ap 2 ) 5 /2 

2 +P ° + lXP ( p 2+Ap2)3/2 



(24) 



is shown in Fig. [2] It is seen that the approximate so- 



long-wavelength (k < 2) part of the spectrum only. Os- 
cillations with shorter wavelength falls in the region of 
strong Landau damping and demonstrate the anomalous 
dispersion (duj r /dk < 0). For ion-acoustic waves, on the 
contrary, the kinetic approximation ( 19 ) turns out to be 



so accurate that both real and imaginary parts of com- 
plex frequency fit the numerical solution of exact disper- 
sion equation in the wide range of wavenumbers (fig. |3| . 
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Figure 3: The dispersion (a) and Landau damping (b) of ion- 
acoustic waves in the non-maxwellian plasma (numerical so- 
lutions of equation ( 22 l coincide with the results of kinetic 
approximation) . 



III. MODULATIONAL INSTABILITY 



W 0,,5 


1 1 


1 1 1 1 1 1 ■ 

exact — w 


s & 0,10 






U 0,05 




^ : 


0,00 




hydro 


1 1 1 1 



2 

kc/co r 



Figure 2: The dispersion (a) and Landau damping (b) of 
Langmuir waves in the non-maxwellian plasma. 



Let us consider stability of the monochromatic Lang- 
muir wave 

E(t, r) = X - (Eoe-^ + lkr + c.c.) (25) 

with respect to the coupled system of low-frequency and 
high-frequency modulational perturbations 

SE{t , r ) - ™ e - int + ^ r + ^± e -*w+* + 

+ S ^ e - iuJ - t + ik - T +c.c, (26) 

where 



lution for the parameters of interest is adequate in the 



uj± = Wfc ± O, k± = k ± q. 
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In the plasma with cold ions and arbitrarily distributed 
electrons, equations for the amplitudes of potential high- 
frequency satellites can be written in the following form: 



e+8E+ 

E* 6E + 
kk+ 

e*_5E*_ 
E* 5E + 



where 



Gf = Aire 2 



e 2 E^ 
Akk + 

G4. 

q 2 £ q 



Akk- 

q 2 £ q 



G 



E SE*_ /G+G, 



kk- 



Gi 



(27) 



E Q SE*_ (G^G 



kk^ 



q 2 z q 



GZ 



(28) 



d 3 p 



k+v 



Gf = Aire* 



fl — qv 

dp — x 

lo± — k±v 



\Uk - kv 



(29) 



<9 P / k • d p f 



k-a, 



Q — qv — kv J fi — qv 

j3 k ' d P 

dp ' — x 

LO± — k±V 

k_ • d p f \ k ■ d n 



G^ = Aire 2 



k-a, 



Gi = A-ne 



k 4. • d„ 



Gk = Aire' 



k-a n 



qv \<jJ— — k_v 
d 3 p 



-k. 



(30) 



— qv ywfe — kv 



SI — qv 

w fc - kv 

d 3 p 
SI — qv 

k- ■ dp/ 

j k_\ 



-k-a 



f k+ ■ d p / 
p \w+ — k + v 



/ k ■ gp/ 

V^fe - kv 



(31) 



(32) 



(33) 



£± = £||(w ± ,k±), £g = £||(n,q). 

In the limiting cases SI <C qv and u>k,± S> k±v, all inte- 
grals Gi reduce to averaging < 1/v 2 >. Thus, the main 
contribution to these integrals comes from low-energy 
core electrons, the typical velocity of which is really much 
smaller than phase velocities of high-frequency satellites 
and the speed of light. It allows us to neglect relativistic 
effects in calculations of G; and transform the equations 
to the well-known form |2"01 



E Q cos ( 



^ (E*5E + cost 



E SE*_ cos 6L 



(34) 



et6El=- E ° COs9 - 



16irnT« 



uljtt 2 , cos9± = 



(E^SE+ cos 9+ + E Q SE1 cos 6_) 

(35) 
(36) 



kk d 

kk-i 



The specihc features of non-maxwellian plasma in these 
equations manifest themselves only in changing of tem- 
perature notation from T to T s and in accurate calcula- 
tion of permittivities e± in the framework of relativistic 
kinetic theory. 

Thus, the modulational instability in a non-maxwellian 
plasma can be described by the following equation: 



n 2 ~n 2 



cos 2 9 + cos 2 



where 



si 



q 2 T s 



(37) 



1 + q 2 T s 



corresponds to the linear dispersion of ion-acoustic waves, 
and the parameter 



W = 



lQimT s 

determines the relative energy of the Langmuir pump 
wave. Since we look for unstable solutions of (37) SI = 
Re O + «Im $7 with Im SI > Im ujk , the last term in the 
definition of £y , appearing as a result of analytical con- 
tinuation of this function to the lower half-plane of w, 
must be omitted. 

For example, let us find a numerical solution of the 



dispersion relation ( 37 1 in the case when the role of long- 
wavelength pump is played by the undamped Langmuir 
wave with the wavenumber k = 1, the frequency w& = 
1.00203 and the energy W = 0.01. Fig. g] shows that 
the growth rate of modulational instability Imf2(gx,gii) 
reaches the maximum for the perturbations propagating 
along k. Since this maximum falls in the region of rather 
small wavenumbers gy ~ 2, it is reasonable to estimate 
how accurately the unstable spectrum can be described 
by the fluid approach, in which thermal corrections to the 
dispersion of high-frequency oscillations depends on the 
increased effective temperature Tg > T s . In this limit, 
the functions e± in the dispersion equation should have 
the simple form 



e± = 1- 



1 



3k 2 ± T e 



(38) 



The unstable spectrum of longitudinally propagating 
modulational perturbations (qj_ = 0) in comparison with 
the case of maxwellian plasma (Tg = T s ) is shown in Fig. 
[5] It is seen that the increase in the temperature Tg with 
respect to T s in the hydrodynamic limit results in almost 
the same narrowing of the spectrum of modulational in- 
stability (by the factor Tg/T s ) that is observed in the case 
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Figure 4: The growth-rate map for the modulational instabil- 
ity pumped by a monochromatic Langmuir wave. 



of exact kinetic description of Langmuir waves. Another 
important result is that the appearance of the tail compo- 
nent in the electron distribution does not change neither 
the maximum of the growth rate nor its position in the 
wavenumber space. It means that we can estimate the 
typical gtowth rate and the corresponding wavenumber 
of modulational instability in the non-maxwellian plasma 
without the detailed study of the real distribution func- 
tion. For this purpose, we can restrict ourselves by the 
case of maxwellian plasma with the temperature corre- 
sponding to the temperature of core electrons. 



IV. SUMMARY 

In the present study we have calculated dispersion laws 
and damping rates of Langmuir and ion-acoustic oscil- 
lations in the case, when an isotropic plasma has the 
energetic tail of superthermal electrons typical to the 
beam-plasma experiments at the multi-mirror trap GOL- 
3. The specific feature of such distribution is that most 
of the kinetic energy is concentrated in a small popula- 
tion of fast electrons. The concept of temperature in this 
case loses its usual meaning, since the energy spread of 
the bulk of electrons appears to be much smaller than 
the total kinetic energy. Nevertheless, the dispersion of 
linear Langmuir waves in the long-wavelength region can 
be correctly described by the fluid approach, in which 
the effective temperature Tn exceeds substantially the 
temperature of core electrons T c , but, due to relativis- 
ts effects, does not reach the temperature of the whole 
distribution. Ion-acoustic waves, in turn, are shown to 
be well described by the kinetic approximation with the 
temperature of low-energy core electrons T s ~ T c . 

Modifications in the linear dispersion of high-frequency 
waves in a non-maxwellian plasma should have an im- 
pact on the build-up efficiency of modulational insta- 
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Figure 5: The real (a) and imaginary (b) parts of the fre- 
quency of longitudinally propagating modulational perturba- 
tions (the fluid approximation for the maxwellian plasma with 
Tt — T s - dashed line, the fluid approximation for the non- 
maxwellian plasma with the increased temperature Te - blue 



solid line, numerical solution of (371 with the exact kinetic 
expressions for e± - red solid line)! 



bility playing the important role in the scenarios of 
strong plasma turbulence. We have solved numerically 
the dispersion equation for modulational perturbations, 
in which both kinetic and relativistic effects are taken 
into account. It is shown that the unstable spectrum in 
the non-maxwellian plasma narrows significantly without 
changing in the growth rate and wavelength of the most 
unstable perturbation. It is also found that the same 
trend is qualitatively reproduced in the case when high- 
frequency waves are described by the fluid approximation 
with the increased effective temperature Ti> T s . 

Independence of the typical growth rate and the wave- 
length of modulational instability on distribution details 
of superthermal electrons allows to prove the possibility 
of applying the model of strong Langmuir turbulence [12] , 
used for calculation of electromagnetic plasma emission, 
to the case of realistic non-maxwellian plasmas typical to 
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laboratory beam-plasma experiments. 
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